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\q 1 I. INTRODUCTION 

O; 

Our universe is known to be homogeneous and isotropic to a very high degree of precision [11,1. Such an universe 
can be described by the well known Friedmann-Robertson- Walker (FRW) metric 0]. There had been, however, 
some cosmological problems associated with the standard big bang model responsible for the evolution of our present 
universe. Inflationary models resolves many problems including the flatness, monopole, and horizon problem^. 
U ' Moreover, gravitational physics is expected to be different from the Einstein-Hilbert models near the Planck scale 
fC) , 0- @|- For example, quantum gravity or string corrections indicate that higher derivative terms could have some 
■ interesting cosmological applicastions [5( . In addition, higher derivative terms also arise as the quantum corrections of 
the matter fields Q. Therefore, the possibility of deriving inflation from higher derivative corrections have been an 
t-H . focus of research interest fH M. M. [Tol. ITiT] . 

For example, a general analysis on the stability conditions of gravity theories has been shown to be useful in the 
search of physical models compatible with our physical universe. In particular, the stability condition for a variety of 
pure gravity theories as a potential candidate of inflationary universe in the flat Friedmann-Robertson- Walker (FRW) 
! space is derived in Ref . [HI, [12, EH . 

— ■ In addition, the highly isotropic universe should be stable against anisotropic perturbation so that the isotropic 
\Q ' FRW space can remain a stable final state. The initial state of our universe could also be anisotropic before it 

becomes isotropic at later stage of the evolution. Nonetheless, it is interesting itself to study the stability analysis 
of the anisotropic space during the post-inflationary epoch. For example, it is shown that anisotropic inflationary 
solution does exist for an NS-NS model with a metric, a dilaton, and an axion field [l4j . This inflationary solution is 

O ' also shown to be stable against small perturbations [lj| . Similar stability analysis has also been studied for a various 

D , of interesting models (l6|. 

Recently, there are also growing interests in the study of Kantowski-Sachs (KS) type anisotropic spaces (Th. fljl . 
We will hence try to study the problem of existence and stability associated with the inflationary solution with an de 
k~j ' Sitter final state. In particular, we will study the effects of the higher derivative terms in Kantowski-Sachs spaces. 
In fact, a large class of pure gravity models and induced gravity models with inflationary KS/FRW solutions was 
presented in Ref. [2(| . Any KS type solution leading itself to an asymptotic FRW final state will be referred to as 

- - ' the KS/FRW solution in this paper for convenience. 

It is shown that the stability of the de Sitter background space is closely related to the choice of the coupling 
constants [13]. Indeed, the pure gravity model given below: 

C g = -R- aR 2 - (3RZR» + 7 i?% 7 R^l p R^ v (I) 

admits an inflationary solution with a constant Hubble parameter given by Hq = l/4r/. This will requires that 7 > 0. 
Here a, /3, and 7 are coupling constants. This shows that (a) the 7 factor determines the scale of the inflation 
characterized by the Hubble parameter Hq and (b) a and (3 factors are irrelevant to the scale Hq in the de Sitter 
phase. The quadratic terms are, however, important to the stability of the de Sitter phase. 
Indeed, perturbing the KS type metric with Hi — ► Hq + SHi, we can show that 

SHi = Cl exp[=^(l + 5,)] + di expp^(l - S,)} (2) 

for 

<*i = y / l + 8/[27-9(6a + 2/3)iJ 2 ] (3) 
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and some arbitrary constants c$, di to be determined by the initial perturbations. Here Hi = dj/aj with aj(i) the scale 
factor in i-direction. We will describe the notation shortly in section II. It is easy to see that any small perturbation 
8Hi will be stable against the de Sitter background if both modes characterized by the exponents 

A± = -[3flot/2][l±«i] (4) 

are all negative. This will happen if S\ < 1. In such case, the inflationary de Sitter space will remain a stable 
background as the universe evolves. 

It can be shown that the stability equation for the anisotropic KS space and the stability equation for the isotropic 
FRW space in the presence of the same inflationary de Sitter background turns out to be identical [ll], EH EH- 
Therefore, the stability of isotropic perturbations also ensures the stability of the anisotropic perturbations. The 
stability of the isotropic perturbations for the FRW space is important for any physical models. Unfortunately, 
inflationary models that are stable against any isotropic perturbations will have problem with the graceful exit 
process. Therefore, the pure gravity model may have troubles dealing with the stability and exit mechanism all 
together. 

Instead of the pure gravity theory, a slow rollover scalar field may help resolving this problem. An inflationary de 
Sitter solution in a scalar-tensor model is expected to have one stable mode (against the perturbation in 8Hi direction) 
and one unstable mode (against the perturbation in 5<j> direction). As a result, the inflationary era will come to an 
end once the unstable mode takes over after a brief period of inflationary expansion. Therefore, we propose to study 
the effect of such theory. 

In particular, we will show in this paper that the roles played by the higher derivative terms are dramatically 
different in the inflationary phase of our physical universe in both pure gravity theory and scalar-tensor theory. First 
of all, third order term will be shown to determine the expansion rate Ho for the inflationary de Sitter space. The 
quadratic terms will be shown to have nothing to do with the expansion rate of the background de Sitter space. They 
will however affect the stability condition of the de Sitter phase. Their roles played in the existence and stability 
condition of the evolution of the de Sitter space are dramatically different. 



II. NON-REDUNDANT FIELD EQUATION AND BIANCHI IDENTITY IN KS SPACE 



Given the metric of the following form: 

ds 2 = -dt 2 + c 2 (t)dr 2 + a 2 (t){d 2 d + f 2 {6)dtp 2 ) (5) 

with f(9) — (9, sinh 9, sin#) denoting the flat, open and close anisotropic space known as Kantowski-Sachs type 
anisotropic spaces. More specifically, Bianchi I (BI), III(BIII), and Kantowski-Sachs (KS) space corresponds to the 
flat, open and closed model respectively. This metric can be rewritten as 

dr 2 

ds 2 = -dt 2 + a 2 {t){ — + r 2 d9 2 ) + a 2 z (t)dz 2 (6) 

1 — kr z 

with r, 9, and z read as the polar coordinates and z coordinate for convenience and for easier comparison with the 
FRW metric. Note that k — 0, 1, —1 stands for the flat, open and closed universes similar to the FRW space. 

Writing H^ = G^ v — T^ Ul Einstein equation can be written as D^H 111 ' — incorporating the Bianchi identity 
D^G^" — and the energy momentum conservation D^T^ = 0. Here G MI/ and represent the Einstein tensor 
and the energy momentum tensor coupled to the system respectively. With the metric ([S]), it can be shown that the 
r component of the equation D^H^ = implies that 

H\ = H%. (7) 

This result also says that any matter coupled to the system has the symmetric property T r r = T e g . In addition, the 
equations D^H^ = and D^H^ Z = both vanish identically for all kinds of energy momentum tensors. More 
interesting information comes from the t component of this equation. It says: 

(d t + 3H)H\ = 2H 1 H r r + H Z H Z Z . (8) 

This equation implies that (i) H\ = implies that H r r = H z z = and (ii) H r r = H z z = only implies (d t + 3H)H t t = 
instead of H\ = 0. Case (ii) can be solved to give H\ = constant x exp[— a 2 a z ] which approaches zero when a 2 a z — > oo. 
For the anisotropic KS spaces, the metric contains two independent variables a and a z . The Einstein field equations 
have, however, three non-vanishing components: H\ = 0, H r r = H e g = and H z z = 0. The Bianchi identity implies 
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that the tt component is not redundant and will hence be retained for complete analysis. Ignoring either one of the 
rr or zz components will not affect the final result of the system. In short, the H\ — equation, known as the 
generalized Friedman equation, is a non-redundant field equation as compared to the H r r = and H* = equations. 

In addition, restoring the g tt component b 2 (t) = 1/B\ will be helpful in deriving the non-redundant field equation 
associated with Gu that will be shown shortly. More specifically, the generalized KS metric will be written as: 

dv 2 

ds 2 = -b 2 (t)dt 2 + a 2 {t){- — + r 2 d6 2 ) + a\[t)dz 2 . (9) 

In principle, the Lagrangian of the system can be reduced from a functional of the metric g^, £((^„), to a simpler 
function of a(t) and a z (t), namely Lit) = a 2 a z L{g^ v {a{t),a z {t))). The equation of motion should be reconstructed 
from the variation of the reduced Lagrangian L(t) with respect to the variable a and a z . The result is, however, 
incomplete because, the variation of a and a z are related to the variation of g rr and g zz respectively. The field 
equation from varying gu can not be derived without restoring the variable b(t) in advance. This is the motivation 
to introduce the metric © such that the reduced Lagrangian L(t) = ba 2 a z C(g^(b(t),a(t),a z (t))) retains the non- 
redundant information of the Hj = equation. Non-redundant Friedmann equation can be reproduced resetting b = 1 
after the variation of b(t) has been done. 

After some algebra, all non-vanishing components of the curvature tensor can be computed: [2p| 

R u tj = ^B 1 H l + B 1 {H t +H 2 )}S} : (10) 



R l \ t = B 1 H i H j e^ m e klm + 3^e il2 (11) 



with Hi = (a/a,a/a,d z /a. z ) = {H\,H% = H\,H Z ) for r, 8, and z component respectively. 
Given a Lagrangian L — ^fgC = L(b(t), a{(t), a z (t)), it can be shown that 



a 2 a z „ , „ a 2 a 



^(#% i^) = ^C{H u H it c?) (12) 



The variational equations for this action can be shown to be: (20| 

C + H t (^- + 3H)L l = HtLt + HiL* (13) 
at 

£ + (| + 3H)»L- = (f t+ 3H)L z (14) 

Here Lj = SC/SHi, L l = 8C/8Hi, and 3H = Hi. For simplicity, we will write C as L from now on in this paper. 
As a result, the field equations can be written in a more comprehensive form: 

DL = L + Hi^ + m)^ -HiU-HiV = (15) 

D Z L = L+{j t +'SH) 2 L z - &+3H)L z = (16) 

III. HIGHER DERIVATIVE GRAVITY MODEL WITH A SCALAR FIELD 

In this section, we will study the higher derivative induced gravity model: 

L = —R - aR 2 — + 7#% 7 R P l P R a % - \d^cp - V{<j>) = L g + L<p (17) 

with L g = —R— aR 2 - f3R^R^ + R""^ R & \ p R ap pv and L^ = -^d^d^cf) - V(<f>) denoting the pure gravity terms and 
the scalar field Lagrangian respectively. 

The corresponding Lagrangian can be shown to be: 

L = 2(2A + B + 2C + D)-Aa [4A 2 + B 2 + AC 2 + D 2 + AAB + 8AC + AAD + ABC + 2BD + 4CT>] 
-2(3 [3A 2 + B 2 + 3C 2 + D 2 + 2AB + 2AC + 2AD + 2BC + 2CD] + 87 [2A 3 + B 3 + 2C 3 + D 3 ] 

+ - V(4>) (18) 
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with 



A = Hi + Hi 


(19) 


B = Hf + ^, 

a z 


(20) 


C = HiH zl 


(21) 


D = H Z + H 2 Z . 


(22) 



This Lagrangian can be shown to reproduce the de Sitter models when we set Hi — > H in the isotropic limit. The 
Friedmann equation reads: 

DL g = \fi + V{<j>) (23) 

for the induced gravity model. In addition, the scalar field equation can be shown to be: 

+ 3Ho4> + V' = 0. (24) 

The leading order de Sitter solution with </> = </> and Hi = H n for all directions can be shown to be: 

V = V(M =6[l-4 7 H*]Hl (25) 

Note that the coupling constants a and (3 do not affect the strength of inflation determined by the Hubble constant 
Hq. This polynomial equation can be written as an cubic equation in terms of the variable x = H^: 

/( ,)^3_i_, + A = o. (26) 

It can be shown that the cubic polynomial f(x) attains its maximum and minimum at x = xm and x = x m respectively 
with x m = —xm = 1/[2 v / 3t]. Plotting C = (x, y = f(x) ) as a curve C on the x-y plane, it is easy to show that the 
curve C intersects with the y-axis at the point Y = (0, Vo/[24 7 ]). 

If f(y) attains its maximum and minimum at the points M — (%,/(im)) an d = (x mi f(x m )) respectively, it 
can be shown that f(x = ±x m ) = (Vo =F ^ / V^l) / [^l] ■ Note that the function f{x) — > ±oo asu ±oo. Therefore 
curve C will have two intersection points x + = (x + , 0) and x_ = (x_, 0) with the positive x-axis if the the minimum 
point N locates at the quadrant IV, or equivalently, f(x m ) < 0. This implies that the cubic equation f(x) = has two 
positive roots if 3"fVQ < 4. There is, however, only one positive root when x+ = x_, or equivalently, when 3jVq = 4. 

Indeed, the cubic equation f(x) — can be solved by using the triple angle formula of cosine function: 

cos[30] =4 cos 3 (9 -3 cos (9 . (27) 
As a result, the solutions to the cubic function can be shown to be 



l 



1/2 



{*;) cos l^-J- (29) 

with defined by the identity cos^o = v^7 Vo/2 < 1. The notation x± is defined according to its orientation with 
respect to the x-coordinate, i.e. x_ < x m < x + . These solutions can be easily converted to the well known solutions 
when 37V 2 > 4 remains valid. 

Note that H = x^}^ 2 > can have two different choices as long as i^V^ < 4 is satisfied. These two solutions 
become identical to each other when 37 V 2 = 4, or equivalently, when cos#o = 1- It is also straightforward to verify 
that Xi given above does solve the cubic polynomial equation f(x) — in consistent with the constraint: 37V 2 < 4. 
In addition, x\ can not be a physical solution because it is negative. 



IV. STABILITY OF HIGHER DERIVATIVE INFLATIONARY SOLUTION 



Our universe could start out anisotropic and evolves to the present highly isotropic state in the post inflationary 
era. Therefore, a stable KS/FRW solution is necessary for any physical model of our universe. 
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Given an effective action of the sort described by Eq. (|45p . the scale factor Hq is determined by the solution 
Hq = x± in the presence of the de Sitter solution Hi = H and the static condition <\> = <f> . In addition to these 
constraints, small perturbations, Hi = H i0 + 5Hi and 4> — 4>q + Scf), against the background de Sitter solution (H i0 , fa) 
may also put a few more constraints to the system. This perturbation will enable us to understand whether the 
background solution is stable or not. 

It can be shown that the perturbation equation of the Bianchi models are identical to the perturbation equation 
of the FRW models. Therefore, any inflationary solutions with a stable mode and an unstable mode could provide 
us with a natural resolution to problem of graceful exit. Such models will, however, also be unstable against the 
anisotropic perturbations. Therefore, any inflationary solutions with a stable mode and an unstable mode is also 
negative to our search for a stable and isotropic inflationary model. It will be shown shortly that the higher derivative 
gravity theory with an inflaton scalar field could hopefully resolve this problem all together. 

In practice, perturbing the background de Sitter solution along the SHi direction should be stable for at least 
a brief moment such that around 60 e-fold inflation can be induced before the de Sitter phase collapses. And the 
resulting universe should also be stable against isotropic and anisotropic perturbations. In addition, the scalar field 
is expected to roll slowly from the initial state cj> = 4>q. Therefore, the perturbation along the 8cj> direction is expected 
to be unstable favoring the system for a natural mechanism of graceful exit. Hence, we will try to study the stability 
equations of the system for small perturbations against the de Sitter background solutions in this section. 

The first order perturbation equation for DL, with Hi — > Ho + SHi, can be shown to be [20l |: 

S(DL) = < H l L ij 8Hj > +3H < > +5 < H t V > +3H < (H.L) + L 3 )SH 3 > 

+ < HiL* > 8(3H)- < HiLijSHj > (30) 

for any DL defined by Eq. <p~5|) with all functions of Hi evaluated at some FRW background Hi = Hq . The notation 
< AiBi >= ^2 i=1 z AiBi denotes the summation over i = 1 and z for repeated indices. Note that we have absorbed 
the information of i = 2 into i = 1 since they contributes equally to the field equations in the KS type spaces. In 
addition, Lj = S 2 L/5HiSHj and similarly for Lij and with upper index % and lower index j denoting variation with 

respect to £f, and Hj respectively for convenience. In addition, perturbing Eq. (|16|) can also be shown to reproduce 
the Eq. (3DJ> in the de Sitter phase 0. 
In addition, it can be shown that 



< H t L a > = 2 < HV Z >, (31) 

< HiL\ > = 2 < HiL\ >, (32) 

L 1 = 2L Z , (33) 

< HiL a > = 2 < H,L lz >, (34) 



in the inflationary de Sitter background with Hq — constant. Therefore, the stability equations (|30|) can be greatly 
simplified. For convenience, we will define the operator T>l as 

V L 5H =< H t L l1 > SH + 3H < H t L l1 > 5H + 3H < H t L\ + L 1 > SH + 2 < H t L l > 5H- < H t L a > 5H. (35) 

As a result, the stability equation ([50)) becomes 

5(DL) = V L {SHi + SH z /2) = ^V L {SH) (36) 

with H = (2Hi + H z )/3 as the average of Hi. 

Hence the leading order perturbation equation in 5H and 8<j) of the Friedmann equation can be shown to be: 

S(DL g ) = ^V g 5H = V'54> (37) 

with T> g SH = T>L g 5H as short-handed notations. This equation can further be shown to be: 

12H 1 2 [67-ffo -3a-0\ {5H + 3H SH) + [l - 12jH^] 5h} = VqSc/). (38) 

Similarly, the leading perturbation of the scalar field equation can be shown to be: 



6$ + 3H a 5cp + V "5cp = 



(39) 
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The variational equation of a z can also be shown explicitly to be redundant in the limit H{ = Hn + SHi and <j) — (f>o + 8(f) 
following the Bianchi identity. 

Assuming that 8H = exp[hHot]8Hn and 8(f> = exp\pHnt]8(f)o for some constants h and p, one can write above 
equations as: 

12H {2 [6 7 # 2 - 3a - 0\ H 2 (h 2 + 3h)8H + [l - 12 7 ff 4 ] 8H } = Vj8<p. (40) 



H 2 (p 2 + 3p+%)80 = 0. (41) 



These equations are consistent when all coefficients vanish simultaneously. This implies that V ' = and 



1 - 127gp 4 
2H 2 [6jH 2 -3a-0\ 



(P 2 + 3P+S) - 0. (43) 



As a result, the solutions to p are p = p± = — 3/2 ± yQ — 4Vj 3 "/_ff 2 /2. Since Vq = 0, the scalar field equation can be 
solved in the de Sitter background Hi = Ho. Indeed, the solution to the equation <j> + 3Ho(f> ~ is: 

"" r l-cxp(-3ff t)]. (44) 



3H, 



This result indicates that the scalar field does change very slowly similar to the slow roll-over assumption in various 
scalar field models. 

An appropriate effective spontaneously symmetry breaking (SSB) potential V of the following form 

V(<t>) = \{tf~v 2 ) 2 + V m (45) 

with arbitrary coupling constant A can be shown to be a good candidate for our physical universe. Here V m is a small 
cosmological constant dressing the SSB potential. The local extremum of this effective potential can be shown to be 
<f> = (local maximum) and <fi = v (local minimum) . We expect that the scalar field starts off from the initial state 
= and rolls slowly down to its local minimum locates at </> = v. 

When the scalar field eventually rolls down to the minimum of V at <f> = v, the system will oscillate around this 
local minimum with a friction term related to the effective Hubble constant H m at this stage. Reheating process 
is expected to take away the kinetic energy of the scalar field. The scalar field will eventually becomes a constant 
background field losing all its dynamics. 

The value of Hq can be chosen to induce enough inflation for a brief moment as long as the slowly changing 
scalar field remains close to the initial state (j> — 4>o — 0. The de Sitter phase will hence remain valid and drive the 
inflationary process for a brief moment determined by the decaying speed of the scalar field. 

In addition, at the initial stage, 



P± = -3/2 ± yj 9 + A\v 2 /H 2 /2 (46) 

for the SSB 4 potential model. Therefore, it is easy to prove that p + > and p_ < indicating an unstable mode 
and a stable mode do exist when we perturb the scalar field in this model. Hence the system will be unstable against 
the perturbation of the scalar field consistent with the slowly changing scalar field solution (j4~4")) . 
On the other hand, the solutions to the /i-equation are h = h± = —3(1 ± 82) /2 with 

81 = 1 + 2(12 7 H 4 - l)/{9ff 2 [6 7 ff 2 - 3a - /?]} (47) 

indicating that h + < 0. Therefore, the /i-perturbation can have two stable modes only when h- < or 8\ < 0. This 
implies that either 

(a) 1 < 12 7 # 4 < 2(3a + (3)H 2 (48) 
or (b) 2(3a + p)H$ < 12 7 i$ < 1 (49) 

has to hold. Let us define the function T± = 12~/H± — 1 with H\ = x± — [cos(6> =p 7r)/3]/y / 3 7 for arbitrary 9. T± 
function can be written as 

r ± = ±4sin-cos[-T|]- (50) 



7 



Therefore, it is easy to show that r + > and T_ < for all < 9 < ir/2. As a result, solution x + is consistent only 
with the constraint (a) in Eq. (|48[) . Similarly, the solution is consistent only with the constraint (b) in Eq. (|4"9"1) . 

In addition, the equation for the perturbation on anisotropic KS type space and isotropic FRW space are the same. 
Therefore, the stable /i-perturbation also ensures that isotropy of the de Sitter space can be made stable against any 
small anisotropic perturbation Hi = Hq + 8Hi . 

Once the scalar field falls close to its local minimum at (f> = v, the scalar field will start to oscillate around this 
local minimum. The Hubble parameter at this stage will read Hi = H rn with H m given by the solution to the leading 
order Friedmann equation at this stage with <j) — v and Hi = H m : 

V m ee V(v) = 6[l-4^]Hl. (51) 

The cubic equation in x = H^ also reads: 

g (x) = x 3 - ^-x +7^ = ( x - x ^)( x ~ x 2+){v - x 2-) = °j ( 52 ) 

with 

1 \ V2 



" 2 = "UJ C "T' (53) 

and 9 m defined by the identity cosOfji = y^3^y V m /2 ^ 1. Note again that only positive roots X2± can be physical 
solutions to H^. In addition, Hq should be much larger than H m in order to generate inflation. Moreover, any 
physical solution should also respect the fact that Vq > V m such that these solutions can be made consistent with the 
SSB potential discussed here. 

Indeed, the function y± = (3 7 ) 1/2 a; ± = cos (9 T tt)/3 (with cos6> = y/lpy V/2) can be shown to be a monotonically 
increasing/decreasing function in 9. Therefore, when 9 decreases, corresponding to the increasing of V, y± will 
decrease/increase accordingly. Hence the relations H$ > H^ and Vq > V m are consistent if we choose the set of 
solutions as (x_, X2-)- Therefore, the Hubble constant H m at the final stage can be chosen to be much smaller than 
the Hubble constant Hq at the inflationary era with proper choice of coupling constants. 

We can also choose the set of solutions as (x+, X2-)- In any case, the initial condition 9q ~ and the final condition 
0m ~ 7r/2 can be adjusted with the proper choice of Vo and V m . Note that #o and 9 m can be chosen to be close to 
and 7r/2 respectively. 

In addition, the first order perturbation equations also have two modes: p — p± — —3/2 ± [9 — AV^/H^ ] 1 / 2 /2 = 
-3/2 ± [9 - 8\v 2 /H? n ] 1 / 2 /2 for the SSB <p 4 potential model gSD. Stability is also ensured even if the discriminator 
becomes imaginary. Therefore, there is no further constraints to be imposed on the parameters from the p-mode 
perturbation. 

There are also two modes for the /i-equation: h = h± = —3(1 ± S 3 )/2 with 

Sj = 1 + 2(12 7 < - 1)/ W 2 n [6 7 i^ - 3a - /?]}. (55) 

Both ft,-modes must be stable for the existence of a stable isotropic de Sitter space in the post-inflationary era. This 
implies that either 

(a) K12 7 ^<2(3a + /3) J ff 2 n (56) 
or (b) 2(3a + /3) J ff 2 n <12 7 <<l (57) 

has to hold. Note that constraint (a) and constraint (b) are organized in a similar way for both Hq and H m . This 
will make our analysis and comparison more easily. 

If we take the set of solutions as (x_,X2-) consistent with Hq ^> H m , both Hq and H m should obey the constraint 
(b). Indeed, Eq.s (|48H49[) imply that 12^/Hq < 1 and \2^Hf n < 1. Similarly, if we take the set of solutions as 
(x + ,X2-), it can be shown that \2^Hq > 1 and 12*yHf n < 1. Therefore, we should choose constraint (a) for Hq 
and constraint (b) for H m . Explicit calculation shows, however, that this will lead to inconsistent constraint on the 
factor ki ee (3a + /3)/(6 7 ). Indeed, the inequality \2^Hq < 2(3a + /3)Hq implies that Hq < k\, but the inequality 
l2-fHf n > 2(3a + (3)H^ implies instead that H? n > k±. This is inconsistent with the fact that Hq > H m . Therefore 
this case can be excluded. As a result, the only consistent set of solutions for Hq and H m are derived from the set of 
solutions (x_,X2_). 
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Consider the special case where V m — 0. The cubic equation becomes x(Ajx 2 — 1) = 0. This implies that 4:jH^ n — 1. 
Or equivalcntly, this solution can also be obtained by taking the limit 8 m = tt/2 in Eq.s ([551 [53}. Therefore, the 
solution to h equation is h = h± = —3(1 ± 5 m )/2 with 

5 m = v / l + 8/[27-9(6a + 2/3)iJ2J. (58) 

The situation at this stage is similar to the case of pure gravity model with a vanishing cosmological constant [20| . 

For the special case V m — 0, we also expect that these two anisotropic perturbation modes against the isotropic 
background solution Hi = H m + 5Hi are both stable modes. As a result, A± = — [3Hot/2][l ± S m ] < can be shown 
to be the necessary condition, implying that S m < 1 or, equivalently, (3a + > 35/18. In particular, both modes 

becomes oscillatory triangular functions when 35/18 > (3a + j3)H^ n > 3/2. In both cases, the de Sitter space will 
remain a stable background as the universe evolves. 



V. CONCLUSION 



The existence of a stable de Sitter background is closely related to the choices of the coupling constants. The 
pure higher derivative gravity model with quadratic and cubic interactions J20| admits an inflationary solution with 
a constant Hubble parameter. Proper choices of the coupling constants allow the de Sitter phase to admit one stable 
mode and one unstable mode for the anisotropic perturbation. 

The stable mode favors a strong inflationary period and the unstable mode provides a natural mechanism for the 
graceful exit process. It is also found that the perturbation against the isotropic FRW background space and the 
perturbation against the anisotropic KS type background space obey the same perturbation equations. This is true 
for both pure and induced gravity models. As a result, the unstable mode in pure gravity model also means that the 
isotropic de Sitter background is unstable against anisotropic perturbations. Therefore, small anisotropic could be 
generated during the de Sitter phase for pure gravity model. 

We have shown that, for gravity models with an additional inflaton scalar field, stable mode for perturbations along 
the anisotropic SHi directions do exist with proper constraints to be imposed on the coupling constants. In addition, 
another unstable mode also exist for the perturbation against the scalar field background </>o with proper constraints. 
Therefore, the de Sitter background can remain stable and isotropic during the inflationary process for this inflaton 
gravity models. 

Explicit model with a spontaneously symmetry breaking (f> potential is presented as an example. Proper constraints 
are derived for reference. It is shown that the these conditions must hold for the initial de Sitter phase solutions in 
order to admit a constant background solution during the inflationary phase for the system. 

Once the scalar field rolls down to the minimum at cf> = v, the system picks up a small Hubble constant H m 
characterized by the small cosmological constant V m . We also show that the resulting de Sitter space can be a stable 
final state with proper choices of coupling constants. Therefore, the scalar-tensor theory proposed in this paper is 
capable of inducing inflation and maintaining a stable FRW de Sitter background space all together. 

In summary, we have shown that a stable mode for (an)isotropic perturbation against the de Sitter background 
does exist for the gravity model with a scalar field. The problem of graceful exit can rely on the unstable mode of 
the scalar field perturbation against the constant phase 4>q. 

It is also found that the quadratic terms will not affect the the inflationary solution characterized by the Hubble 
parameter Hq and H m . These quadratic terms play, however, critical role in the stability of the de Sitter background. 
In addition, it is also interesting to find that their coupling constants a and [3 always show up as a linear combination 
of 3a + (3 in these stability equations. Implications of these constraints deserve more attention for the applications to 
the inflationary models. 
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